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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 


(metric)  units  as  follows: 

Multiply 

By  . 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

kips  (force) 

4.448222 

kilonewtons 
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USER'S  GUIDE:  UTEXAS2  SLOPE-STABILITY  PACKAGE 


PART  I :  INTRODUCTION 


Purpose  of  UTEXAS2 

1.  A  large  number  of  slope-stability  computer  programs  have  been  devel¬ 
oped  over  the  past  years.  A  survey  of  geotechnical  computer  programs  compiled 
by  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  in  1982  (Edris  and 
Vanadit-Ellis  1982)  listed  37  different  programs  used  throughout  the  Corps. 
Being  developed  for  specific  purposes,  many  of  the  programs  are  restricted  in 
the  range  of  conditions  which  can  be  analyzed.  Many  of  the  programs  are  not 
well  documented,  and  the  specific  algorithms  implemented  are  not  readily 
apparent.  Systematic  evaluation  of  the  various  programs  is  difficult  because 
of  their  diversity  and  because  each  requires  a  different  format  for  data 
entry. 

2.  For  these  reasons,  a  joint  venture  task  group  of  the  Computer  Appli¬ 
cations  in  Geotechnical  Engineering  (CAGE)  and  the  Geotechnical  Aspects  of  the 
Computer-Aided  Structural  Engineering  (G-CASE)  projects  was  tasked  by  the 
Headquarters,  US  Army  Corps  of  Engineers  (USACE) ,  to  develop  a  slope-stability 
package  suitable  for  Corps-wide  use.  This  package  of  slope-stability  programs 
will  offer  the  following  benefits  to  the  Corps  of  Engineers  by: 

a.  Providing  documented  material  to  the  design  engineer  as  a  main¬ 
tained  and  updated  part  of  the  Corps'  computer  library. 

b.  Facilitating  review  of  District  work  by  Divisions  and 
architectural-engineering  contract  work  by  both  Corps  Districts 
and  Divisions. 

£.  Enabling  different  analysis  procedures  to  be  conveniently  used 
from  a  common  input  data  file. 

3.  The  criteria  for  this  limit-equilibrium  slope-stability  package  are 
contained  in  Miscellaneous  Paper  GL-85-8  (CAGE  Task  Group  on  Slope  Stability 
1985)  which  concluded  that: 

a.  No  program  in  existence  meets  all  the  criteria  outlined  in  the 
report. 

b.  The  program  UTEXAS  (University  of  Texas  Analysis  of  Slopes) , 
developed  by  Wright  and  Roecker  (1984)  for  the  Texas  Highway 
Department,  most  nearly  meets  all  the  criteria. 
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c.  Capability  and  criteria  modifications  to  UTEXAS  would  be  faster 
and  more  cost  effective  than  to  write  a  new  program. 

The  version  of  UTEXAS  containing  the  additional  capabilities  is  called 
UTEXAS 2. 

Organization  of  User's  Manual 

4.  The  UTEXAS2  User’s  Guide  is  organized  into  three  volumes  to  avoid  a 
large,  cumbersome  report.  The  use  of  separate  volumes  also  provides  for  the 
timely  publication  of  the  user-required  guidelines  of  the  program.  Volume  I 
contains  the  user  guidelines  including  instructions  for  input  and  output, 
illustrative  examples,  search  procedure  recommendations,  and  error  message 
explanations.  Volume  II  of  the  series  contains  the  theory  and  derivations  of 
the  equations  used  in  the  program.  Volume  III  consists  of  problems  illus¬ 
trating  coding  procedures  for  generic  problem  types  and  demonstrating  the 
capabilities  and  versatility  of  the  program. 

Mechanics  of  Procedures 

5.  The  mechanics  of  limit-equilibrium  procedures  are  covered  in  this 
study  of  two-dimensional  slope-stability  analysis.  The  three  force  equilib¬ 
rium  procedures  discussed  in  this  report  include  the  wedge  method  described  in 
Engineer  Manual  (EM)  1110-2-1902  (Headquarters,  Department  of  the  Army  1970) 
and  two  procedures  that  use  the  method  of  slices.  These  force  equilibrium 
procedures  use  the  Corps  of  Engineers  Modified  Swedish  side-force  assumption 
of  parallel  side  forces  at  a  user-specified  inclination,  EM  1110-2-1902 
(Headquarters,  Department  of  the  Army  1970),  or  Lowe  and  Karafiath’s  side- 
force  assumption  (1960).  The  mechanics  of  Bishop's  Simplified  procedure 
(Bishop  1955),  and  Spencer’s  complete  equilibrium  procedure  (Spencer  1967)  are 
also  discussed  as  are  sources  of  potential  error  for  all  methods.  The  various 
cases  of  slope-loading  conditions  (i.e.,  steady  seepage,  sudden  drawdown, 
etc.),  characterization  of  material  properties,  and  internal  water-pressure 
determinations  are  not  covered.  These  topics,  discussed  in  EM  1110-2-1902 
(Headquarters,  Department  of  the  Army  1970),  are  very  important  but  are  beyond 
the  scope  of  this  user's  guide. 
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PART  II:  THEORY  OF  LIMIT-EQUILIBRIUM  PROCEDURES 

Definition  of  the  Factor  of  Safety 

6.  The  f"ctor  of  safety  F*  in  this  program  is  defined  with  respect  to 
shear  stren^r\  as 

F  =  -  (1) 

T 

where 

s  =  available  shear  strength 

t  =  shear  stress  required  for  just-stable,  static  equilibrium 
The  shear  strength  (s)  in  Equation  1  is  expressed  in  terms  of  the  Mohr-Coulomb 
failure  criteria.  In  the  case  of  effective  stress  analyses,  the  shear 
strength  is  expressed  by 

s  =  c  +  (o  -  u)  tan  $  (2) 

where 

c,$  =  shear  strength  parameters  expressed  in  terms  of  effective 
stress 

a  =  total  normal  stress 
u  =  pore  water  pressure 
(a  -  u)  *  effective  normal  stress 

For  total  stress  analyses,  the  shear  strength  parameters  are  expressed  in 
terms  of  total  stress  equivalents,  (c,  <j>)  ,  and  pore  water  pressures  are  not 
considered.  Thus,  total  stress  analyses  are  expressed  by  the  equation 

s  =  c  +  o  tan  ip  (3) 

The  only  differences  between  total  and  effective  stress  expressions  for  shear 
strength  and  factor  of  safety  are  whether  or  not  total  stress  or  effective 
stress-strength  parameters  (c,  $  or  c,  $)  are  used  and  whether  the  pore  water 
pressure  u  appears  in  the  strength  equations.  Procedures  and  equations 
developed  on  the  basis  of  effective  stresses  may  be  applied  equally  to  total 
stress  analyses  by  using  c  and  <\>  ,  rather  than  c  and  $  ,  and  by  con¬ 
sidering  the  pore  water  pressure  to  be  zero. 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  B) . 
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7.  The  factor  of  safety  defined  by  Equation  1  is  computed  from  require¬ 
ments  of  static  equilibrium.  The  factor  of  safety  is  introduced  into  the 
equations  of  static  equilibrium  by  first  expressing  Equation  1  as 


where  t  represents  the  shear  stress  which  appears  in  the  equations  of  static 
equilibrium.  Combining  Equation  4  with  the  expression  for  the  shear  strength 
(Equation  2) ,  the  shear  stress  is  expressed  as 

c  .  ,  .  tan  i  ,r\ 

t  -  -  +  (a  -  u)  — (5) 

where  the  expression  on  the  right-hand  side  of  Equation  5  is  used  to  replace 
the  shear  stress  in  the  equations  of  static  equilibrium.  The  factor  of  safety 
applicable  to  cohesion  (c  ,  c)  is  assumed  to  be  the  same  factor  applicable  to 
the  frictional  component  (tan  $  ,  tan  4>)  of  shear  strength  in  Equation  5.  A 
second  assumption  is  that  the  factor  of  safety  is  constant  along  an  assumed 
shear  (or  sliding)  surface.  These  assumptions  reduce  the  number  of  unknowns 
related  to  the  factor  of  safety  which  must  be  computed  from  the  equilibrium 
equations  and  are  discussed  in  detail  later  in  the  User's  Guide. 

8.  In  limit-equilibrium  procedures,  the  factor  of  safety  is  calculated 
for  an  assumed  shear  surface.  A  number  of  trial  shear  surfaces  are  analyzed 
until  the  one  producing  a  minimum  factory  of  safety  is  found.  In  comparing 
the  factors  of  safety  calculated  by  the  various  procedures,  it  is  appropriate 
to  compare  only  the  minimum  factors.  The  corresponding  critical  shear  sur¬ 
faces  may  differ  for  different  procedures  and,  thus,  separate  critical  shear 
surfaces  must  be  found  for  each  procedure  (Duncan  and  Wright  1980) .  If  fac¬ 
tors  of  safety  for  other  than  the  critical  shear  surfaces  are  compared,  the 
differences  in  results  for  various  procedures  may  be  much  larger  than  the 
differences  among  minimum  factors  of  safety.  The  larger  differences  associ¬ 
ated  with  other-than-critical  shear  surfaces  can  be  misleading  and,  in  most 
cases,  have  no  practical  meaning. 


Subdivision  of  Soil  Mass 


9.  In  order  to  facilitate  limit-equilibrium  procedures,  the  soil  mass 
which  is  bounded  by  the  assumed  shear  surface  and  slope  surface  is  subdivided 
into  a  finite  number  of  vertical  slices.  This  allows  analysis  of  relatively 
inhomogeneous  slopes  where  the  properties  along  the  shear  surface  at  the  base 
of  each  slice  may  differ  from  slice  to  slice. 

10.  Total  forces  acting  on  the  free-body  diagram  of  a  typical  slice 
employed  in  the  procedures  of  slices  are  shown  in  Figure  1.  The  total  forces 


Figure  1.  Total  forces  acting  on  a 
typical  slice  for  the  procedure  of 
slices 


acting  on  the  slice  include  the  weight  of  the  slice  W  ,  the  shear  and  normal 
force  on  the  base  of  the  slice,  S  and  N  ,  respectively,  the  shear  forces  on 
the  left  and  right  boundaries  of  the  slice,  Xi  and  X^,+1  ,  respectively,  and 
the  horizontal  forces  on  the  left  and  right  boundaries  of  the  slice,  and 

E  ,  respectively.  In  the  case  of  effective  stress  analysis,  N  ,  X  ,  and 
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E  may  be  expressed  as  N  ,  X  ,  E  ,  and  a  pore  pressure  term.  In  addition, 
forces  such  as  water  forces  on  the  base  and  sides  ,  and  HR  ;  exter¬ 

nal  surface  loads  P  ;  external  shear  forces  T  ;  and  "pseudo-static"  seismic 
forces  K  •  W  may  also  be  included.  These  forces  are  considered  as  known 
quantities  in  terms  of  their  location,  direction,  and  magnitude.  The  shear 
force  S  on  the  base  of  the  slice  is  defined  in  terms  of  the  normal  force  N 
and  factor  of  safety  F  through  the  Mohr-Coulomb  Strength  Equation  5.  Thus, 
the  shear  force  can  also  be  considered  a  known  quantity  based  on  given  soil 
parameters,  provided  that  F  and  N  are  known  or  calculated  from  equilibrium 
requirements.  The  unknown  quantities  which  must  be  assumed  or  computed  in  a 
solution  of  the  equations  of  static  equilibrium  include  the  F  ,  N  ,  and  the 
side  forces  X  and  E  .  In  addition,  the  locations  of  N  on  the  base  of  the 

slice  JL.  and  of  the  horizontal  side  forces  on  the  sides  of  the  slice  h 
N  t 

are  unknown  and  must  be  assumed  or  computed.  The  total  number  of  unknowns 
associated  with  a  soil  mass  subdivided  into  n  slices  is  summarized  in 
Table  1  (Wright  1982)  with  the  corresponding  number  of  static-equilibrium 
equations;  5n  -  2  unknowns  and  3n  equilibrium  equations  exist. 

Table  1 

Unknowns  and  Equilibrium  Equations  for  Procedures  of  Slices 

(Wright  1982) 


_ Description _  Number* 

Unknowns 

Normal  force  on  base  of  slice  N  n 

Normal  (horizontal)  force  between  slices  E  n  -  1 

Shear  (vertical)  force  between  slices  X  n  -  1 

Location  of  normal  force  on  base  of  slice  n 

N 

Location  of  side  force  between  slices  h  n  -  1 

Factor  of  safety  F  1 

Total  unknowns  5n  -  2 

Equilibrium  equations 

Summation  of  forces  in  the  vertical  direction  n 

Summation  of  forces  in  the  horizontal  direction  n 

Summation  of  moments  n 

Total  equations  3n 


*  Number  of  slices  =  n. 
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11.  The  number  of  unknown  quantities  (5n  -  2)  shown  in  Table  1  exceeds 
the  corresponding  number  of  equilibrium  equations  (3n) .  Accordingly,  assump¬ 
tions  must  be  introduced  to  obtain  a  statically  determinate  solution  for  the 
F  and  for  any  of  the  other  unknown  quantities  which  are  to  be  calculated. 
Various  assumptions  pertaining  to  the  unknown  quantities  are  made  in  each  of 
the  different  limit-equilibrium  procedures.  Differences  among  the  assumptions 
constitute  one  of  the  principal  differences  among  the  various  methods  of 
slices.  In  addition,  some  procedures  satisfy  all  requirements  (3n)  of  planar 
equilibrium  while  others  only  partially  satisfy  complete  equilibrium.  Differ¬ 
ences  among  the  specific  equilibrium  conditions  satisfied  constitute  a  second 
principal  difference  among  the  various  methods  of  slices  procedures. 

Simplified  Bishop  Procedure 

12.  The  Simplified  Bishop,  also  known  as  Bishop's  Modified  procedure 
(Bishop  1955),  is  one  of  the  most  widely  used  and  accepted  methods  of  slices. 
This  procedure,  for  circular  shear  surfaces,  was  proposed  by  Bishop  (1955)  as 
a  simplified  version  of  his  detailed  approach  satisfying  complete  equilibrium. 
According  to  this  procedure,  there  are  assumed  to  be  no  shear  forces  between 
slices  (X  =  0) ,  and  forces  are  resolved  in  the  vertical  direction  to  obtain  an 
equation  for  N  on  the  base  of  each  slice.  The  equation  for  N  is  incor¬ 
porated  into  an  equation  for  moment  equilibrium  about  the  center  of  a  circular 
shear  surface.  The  equation  of  moment  equilibrium  considers  the  entire  soil 
mass  (all  slices)  as  a  single  free  body  and  is  used  to  compute  the  factor  of 
safety.  For  a  slope  with  n  slices,  the  Simplified  Bishop  procedure  employs 
and  satisfies  n  +  1  equilibrium  equations,  consisting  of  n  equations  for 
equilibrium  of  forces  in  the  vertical  direction  and  one  equation  for  moment 
equilibrium  for  the  entire  soil  mass.  The  corresponding  unknowns  which  are 
solved  from  these  equations  are  N  on  the  base  of  each  slice  and  one  factor 
of  safety.  Table  2  summarizes  the  assumptions  along  with  the  unknown  forces 
for  the  Simplified  Bishop  procedure  (Bishop  1955). 

13.  Although  this  procedure  does  not  satisfy  all  requirements  for  sta¬ 
tic  equilibrium,  it  has  been  shown  to  produce  reasonably  correct  values  for 
the  factor  of  safety  (Duncan  and  Wright  1980).  A  number  of  comparisons  for 
relatively  homogeneous  slope  conditions  have  been  made  between  factors  of 
safety  calculated  by  the  Simplified  Bishop  procedure  and  by  procedures  of 
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Assumptions 

Vertical  shear  force  on  sides  of  slice  X 
equals  0 

Horizontal  side  forces  are  not  considered 
Unknowns 

Normal  force  on  base  of  slice  N  n 

Factor  of  safety  F  1 

Total  unknowns  n  +  1 

Equilibrium  equations 

Summation  of  forces  in  the  vertical  direction  n 

Summation  of  moment  of  forces  overall  1 

equilibrium  equation  for  total  soil  mass 

Total  equations  n  +  1 


*  Number  of  slices  =  n. 

slices  which  satisfy  complete  static  equilibrium  (Whitman  and  Bailey  1967, 
Fredlund  and  Krahn  1977,  Duncan  and  Wright  1980).  The  factors  of  safety  are 
usually  in  agreement  within  a  few  percent.  The  factors  calculated  by  the  Sim¬ 
plified  Bishop  procedure  have  also  been  compared  with  values  which  were  calcu¬ 
lated  using  stresses  computed  by  both  linear  and  nonlinear  finite  element 
procedures  (Wright,  Kulhawy,  and  Duncan  1973).  These  comparisons  also  show 
that  the  factors  of  safety  are  in  close  agreement,  ±5  percent  for  homogeneous 
slope  conditions. 

14.  The  principal  limitation  of  the  Simplified  Bishop  procedure  results 
from  the  assumption  of  circular  shear  surfaces.  Although  several  extensions 
of  Bishop's  procedure  to  noncircular  shear  surfaces  have  been  suggested 
(Nonveiller  1965,  Bell  1969),  documentation  of  the  accuracy  of  such  procedures 
is  not  available.  Thus,  the  procedure  is  limited  to  analyses  of  slopes  where 
the  assumption  of  circular  shear  surfaces  is  reasonable.  This  procedure  does 
not  satisfy  horizontal  force  equilibrium,  and  caution  should  be  exercised  when 
used  for  pseudo-static  analyses  of  earthquake  conditions. 


15.  The  derivations  of  the  equations  used  for  the  Simplified  Bishop 
procedure  are  presented  in  Appendix  A. 

Force-Equilibrium  Procedure 

16.  The  force-equilibrium  procedures  are  one  type  of  limit-equilibrium 
procedure  which  consider  and  satisfy  the  requirements  for  equilibrium  of 
forces  in  the  vertical  and  horizontal  direction  for  each  slice  and  the  sliding 
mass,  but  do  not  satisfy  moment  equilibrium.  Such  force-equilibrium  proce¬ 
dures  are  all  suitable  for  analyses  employing  shear  surfaces  of  any  general 
shape,  and  the  required  calculations  can  be  performed  by  hand.  Accordingly, 
force-equilibrium  procedures  have  been  widely  accepted  and  used  by  many  prac¬ 
ticing  engineers.  The  Corps  of  Engineers  uses  this  procedure  as  described  in 
EM  1110-2-1902,  Engineering  and  Design  Stability  of  Earth  and  Rock-Fill  Dams 
(Headquarters,  Department  of  the  Army  1970). 

17.  The  unknowns  for  force-equilibrium  procedures  are  reduced  from  a 

total  of  5n  -  2  to  3n  -  1  by  not  considering  the  locations  £  and  h 

n  t 

associated  with  the  forces  N  and  E  .  There  are  2n  equations  for  force 
equilibrium,  and,  thus,  a  total  of  n  -  1  assumptions  are  required  to  make 
the  problem  statically  determinate.  All  force-equilibrium  procedures  make 
these  n  -  1  assumptions  pertaining  to  the  inclination  of  the  side  forces. 
Side  forces,  X  and  E  ,  can  be  described  in  terms  of  components  of  the 
resultant  Z  and  its  inclination  0  from  the  horizontal.  Thus,  the  side 
forces  shown  in  Figure  1  can  also  be  expressed  by  n  -  1  unknown  values  of  Z 
and  n  -  1  unknown  values  of  0  .  Force-equilibrium  procedures  assume  values 
for  0  to  achieve  a  statically  determinate  solution,  and  Z  is  considered  as 
an  unknown.  The  n  values  for  N  ,  n  -  1  values  for  Z  ,  and  one  value  for 
F  are  computed  by  using  the  2n  equations  of  static  equilibrium  for  all 
slices.  Differences  among  the  many  force-equilibrium  procedures  are  a  result 
of  the  various  assumptions  which  are  made  concerning  the  side-force  inclina¬ 
tion  0  .  As  a  result,  the  factor  of  safety  is  directly  related  to  0  as 
shown  in  Figure  2.*  Table  3  summarizes  the  unknowns  and  the  assumptions  asso¬ 
ciated  with  this  procedure. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 
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FACTOR  OF  SAFETY 


-10  0  +10 

INCLINATION  IN  SIDE  FORCES,  DEGREES 


Figure  2.  Influence  of  the  parallel  side-force  inclination  on  the 
values  of  F  calculated  by  force-equilibrium  and  moment- 
equilibrium  procedures  (after  Wright  1969) 
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Table  3 

Assumptions,  Unknowns,  and  Equilibrium  Equations  for  Force- 
Equilibrium  Procedures 

_ Description _ 

Assumptions 

Inclination  of  resultant  interslice  force  0 
is  specified  and  constant 

Unknowns 

Normal  force  on  base  of  slice  N 

Interslice  normal  force  E  or  resultant  inter¬ 
slice  force  Z 

Factor  of  safety  F 
Total  unknowns 
Equilibrium  equations 

Summation  of  forces  in  the  vertical  direction 
Summation  of  forces  in  the  horizontal  direction 

Total  equations 


*  Number  of  slices  *  n. 

18.  Force-equilibrium  procedures  require  an  assumption  for  0  .  There 
are  several  0  assumptions  that  could  be  used.  EM  1110-2-1902,  Engineering 
and  Design  Stability  of  Earth  and  Rock-Fill  Dams  (Headquarters,  Department  of 
the  Army  1970),  specifies  that  the  side-force  inclination  should  be  parallel 
to  the  average  outer  slope  of  the  embankment  as  shown  in  Figure  2.  This 
assumption  has  been  found  to  produce  factors  of  safety  which  are  as  much  as 
14  percent  higher  than  values  calculated  by  procedures  which  satisfy  complete 
equilibrium  (Wright  1969).  Thus,  the  assumption  is  not  conservative  and  may 
lead  to  significant  overestimates  of  the  factor  of  safety  as  shown  in 
Figure  2. 

19.  There  are  other  0  assumptions  that  could  be  utilized,  the  first 
being  that  side  forces  are  horizontal.  This  assumption  is  the  same  as  that 
employed  in  the  Simplified  Bishop  procedure  (Bishop  1955).  However,  the 
equilibrium  conditions  which  are  satisfied  in  the  force-equilibrium  procedures 
are  different  from  those  satisfied  in  the  Simplified  Bishop  procedure.  The 


Number* 


n 

n  -  1 

1 

2n 


n 

n 

2n 
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assumption  of  horizontal  side  forces  in  the  force-equilibrium  procedures  often 
leads  to  factors  of  safety  which  are  significantly  less  than  those  calculated 
by  the  Simplified  Bishop  procedure  or  procedures  which  satisfy  all  require¬ 
ments  for  static  equilibrium.  Analyses  of  a  variety  of  both  homogeneous  and 
inhomogeneous  slopes  show  that  the  assumption  of  horizontal  side  forces  in 
force-equilibrium  procedures  can  cause  the  factor  of  safety  to  be  at  least 
20  percent  less  than  values  calculated  by  complete  equilibrium  procedures 
(Wright  1969').  Janbu,  Bjerrum,  and  Kjaernsli  (1956)  found  similar  differences 
of  up  to  13  percent. 

20.  Another  side-force  inclination  assumption  uses  the  average  of  the 
outer  slope  and  the  shear-surface  inclination  (Lowe  and  Karafiath  1960).  The 
authors  suggested  that  the  side  force  on  each  boundary  between  slices  acts  at 
an  inclination  which  is  the  average  of  the  inclinations  of  the  slope  and  the 
shear  surface  directly  above  and  below  the  slice  boundary,  respectively.  This 
side-force  assumption  has  been  found  to  produce  values  for  the  factor  of 
safety  which  are  generally  within  ±10  percent  of  the  values  calculated  by 
complete  equilibrium  procedures  (Wright  1969,  Duncan  and  Wright  1980).  The 
differences  (±10  percent)  are  smaller  than  those  associated  with  the  other 
side-force  assumptions  discussed  above.  Thus,  the  user  needs  to  be  aware  of 
the  effect  of  0  assumption  input  for  the  Modified  Swedish  procedure. 

21.  Using  a  0  assumption,  the  Modified  Swedish  procedure  calculates  a 
composite  force  polygon.  The  error  of  closure  is  determined  for  each  trial 
factor  of  safety  until  the  sign  of  the  closure  error  changes.  Then,  the  trial 
factors  of  safety  are  plotted  against  the  error  of  closure,  and  the  equilib¬ 
rium  factor  of  safety  at  zero  error  of  closure  is  determined  by  interpolation. 
Forces  on  typical  slices  and  the  composite-force  polygon  are  shown  in  Fig¬ 
ure  3.  The  derivations  of  the  equations  used  for  the  Modified  Swedish  proce¬ 
dure  are  presented  in  Appendix  A. 

Corps  of  Engineers  Wedge  Procedure 

22.  The  ability  to  analyze  noncircular  shear  surfaces  is  an  advantage 
of  the  force-equilibrium  procedures.  The  Corps  of  Engineers  wedge  procedure 
discussed  in  EM  1110-2-1902  (Headquarters,  Department  of  the  Army  1970)  is  a 
force-equilibrium  procedure  with  specified  side-force  assumptions  that  pro¬ 
vides  the  ability  to  analyze  shear  surfaces  which  correspond  to  zones  of 
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Figure  3.  Modified  Swedish  circular-arc  procedure 


material  weaknesses.  Composite  interslice  or  interblock  forces  are  determined 
for  active  and  passive  wedges  and  central  block.  Typical  force  polygons  are 
shown  in  Figure  4,  and  the  prescribed  inclinations  of  side  forces  are  shown  in 
Figure  5.  These  side-force  inclinations,  where  the  active  and  passive  forces 
are  incited  at  different  angles,  cannot  be  modeled  in  the  program  UTEXAS2. 

23.  The  unknowns,  equilibrium  equations,  and  assumptions  presented  in 
Table  3  for  the  force-equilibrium  procedures  also  apply  to  the  wedge  proce¬ 
dure.  The  derivations  of  the  equations  used  for  the  wedge  procedure  are  the 
same  as  in  Appendix  A  except  for  the  side-force  assumptions  shown  in  Figures  4 
and  5 . 

Limitations  of  the  Force-Equilibrium  Procedures 

24.  The  main  limitation  of  the  force-equilibrium  procedures  results 
from  the  sensitivity  of  the  computed  factor  of  safety  to  the  various  side- 
force  assumptions  since  moment  equilibrium  is  not  considered.  As  shown  in 
Figure  2,  the  factor  of  safety  is  a  function  of  the  side-force  inclination. 
Consequently,  there  is  a  large  degree  of  uncertainty  associated  with  the  cor¬ 
rectness  of  a  given  side-force  assumption  for  a  given  problem.  For  example, 
when  <f>  is  equal  to  zero  and  a  circular  shear  surface  is  assumed,  the  assump¬ 
tion  that  side  forces  are  parallel  to  the  shear  surface  produces  precisely  the 
same  factor  of  safety  as  the  value  which  is  obtained  by  complete  equilibrium 
procedures.  Other  side-force  assumptions,  in  this  case  ( 4>  =0),  produce 
incorrect  results.  However,  when  <p  is  not  equal  to  zero,  the  assumption 
that  side  forces  are  parallel  to  the  shear  surface  may  be  relatively  poor. 

Spencer's  Procedure 

25.  Spencer's  procedure  (Spencer  1967)  is  a  complete  equilibrium  proce¬ 
dure  and  is  defined  as  a  stability-analysis  procedure  that  fully  satisfies  the 
force  and  moment  requirements  of  static  equilibrium  for  each  slice.  Factors 
of  safety  calculated  by  several  of  the  available  complete  equilibrium  proce¬ 
dures  and  compared  for  a  variety  of  different  slopes  (Wright  1969,  Sarma  1973, 
Fredlund  and  Krahn  1977,  Duncan  and  Wright  1980)  have  shown  that  the  results 
of  these  comparisons  are  in  agreement  within  ±5  percent  of  a  mean  value,  and 
no  procedure  appears  to  produce  consistently  high  or  low  values  for  the  factor 
of  safety  (Wright  1982) .  This  procedure  is  suitable  for  analysis  of  both 
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Corps  of  Engineers  wedge  procedure 


Corps  of  Engineers  wedge  method,  prescribed  inclination  of  side  forces 


circular  and  noncircular  or  wedge  shear  surfaces.  However,  the  computations 
require  the  use  of  a  computer  and  are  impractical  for  hand  calculations. 

26.  Spencer  (1967)  developed  this  procedure  which  satisfies  both  force 
and  moment  equilibrium.  The  9  is  not  assumed,  but  determined,  so  that  com¬ 
plete  equilibrium  is  satisfied.  The  side  forces  are  assumed  to  be  parallel. 
Also,  the  assumption  is  made  that  the  normal  forces  are  located  at  the  center 
of  each  slice  base.  These  assumptions  result  in  the  reduction  by  2n  -  2  of 
the  unknowns,  leaving  3n  unknowns.  Thus,  the  problem  becomes  statically 
determinate.  Table  4  summarizes  the  unknowns  and  equilibrium  equations  for 
this  procedure.  The  derivation  of  the  Spencer's  procedure  is  given  in 
Appendix  A. 

27.  It  is  impractical  to  solve  the  equilibrium  equation  in  Spencer's 
procedure  by  hand.  However,  one  can  check  the  method  by  drawing  a  composite- 
force  polygon  in  the  manner  of  the  force-equilibrium  procedure  with  the  Corps 
of  Engineers  Modified  Swedish  side-force  assumption  (EM  1110-2-1902)  (Head¬ 
quarters,  Department  of  the  Army  1970),  using  the  angle  of  side-force  incli¬ 
nation  computed  by  Spencer's  procedure  (Spencer  1967).  The  force  polygon  must 
close  if  the  solution  is  correct.  The  location,  magnitude,  and  sign  of  the 
side  force  for  each  slice  can  be  checked  for  correctness.  Also,  comparison  of 
results  can  be  made  with  the  values  obtained  by  using  one  of  the  other 
procedures . 
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Table  4 

Unknowns,  Assumptions,  and  Equilibrium  Equations  for  the  Complete 
Equilibrium  Procedure  (Spencer  1967) 

_ Description _  Number* 

Assumptions 

Side  forces  assumed  parallel:  X  =  E  •  tan  0  ; 

0  is  the  same  for  all  slices 

Location  of  normal  forces  on  base  of  slice 
assumed  at  midpoint  on  base 

Unknowns 


Normal  force  on  base  of  slice  N  n 

Normal  (horizontal)  force  between  slices  E  n  -  1 

Side  force  inclination  0  1 

Location  of  side  force  between  slices  h  n  -  1 

Factor  of  safety  F  1 

Total  unknowns  3n 

Equilibrium  euations 

Summation  of  forces  in  the  vertical  direction  n 

Summation  of  forces  in  the  horizontal  direction  n 

Summation  of  moments  n 

Total  equations  3n 


*  Number  of  slices  =  n. 
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PART  III:  SOURCES  OF  POTENTIAL  ERROR  (WRIGHT  1982) 


28.  Two  cases  can  be  readily  identified  for  any  of  the  stability  proce¬ 
dures  where  even  the  most  statically  correct  limit-equilibrium  slope-analysis 
procedures  will  produce  incorrect  or  unreasonable  values  for  the  factor  of 
safety  and  for  the  various  other  calculated  values.  The  two  cases  where  in¬ 
correct  or  unreasonable  values  occur  are  related  directly  to  specific  slope 
and  shear  surface  conditions.  The  first  case  where  incorrect  or  unreasonable 
results  may  occur  corresponds  to  slopes  where  material  with  a  relatively  high 
cohesive  component  (c,  c)  of  strength  exists  along  the  upper  crest  portion  of 
the  shear  surface.  The  second  case  corresponds  to  slopes  where  shear  surfaces 
exit  steeply  upward  through  cohesionless  material  near  the  toe  of  the  slope. 
These  two  cases  and  the  resultant  problems  they  introduce  require  separate 
consideration,  and  different  solutions  are  used  to  eliminate  errors  associated 
with  each. 


High  Cohesive-Strength  Components 


29.  When  a  relatively  high  cohesive  component  of  strength  exists  along 
the  shear  surface  near  the  crest  of  a  slope,  limit-equilibrium  solutions  often 
indicate  relatively  high  tensile  stresses;  in  such  cases,  N  acting  on  the 
shear  surface  and  E  near  the  crest  of  the  slope  become  negative  (tensile). 
Also,  near  the  crest  of  the  slope,  the  locations  of  the  side  forces  ("line  of 
thrust")  diverge  and  fall  well  outside  the  boundaries  of  the  slope  and  shear 
surface,  approaching  an  infinite  distance  from  the  slope  in  the  extreme. 

Wright  (1975)  has  shown  that  such  a  pattern  for  the  side-force  locations  is 
closely  related  to  existing  zones  of  tensile  stress.  The  tensile  stresses 
implied  in  the  observed  solutions  are  fully  consistent  with  the  use  of  a  cohe¬ 
sive  component  of  strength  in  the  Mohr-Coulomb  equation  and  the  associated 
stress  state  required  for  a  limit  state  of  equilibrium.  The  zone  in  which  the 
tensile  stresses  are  observed  is  analogous  to  an  active  earth-pressure  zone, 
and  such  tensile  stresses  are  well  recognized  in  classical  active  earth- 
pressure  calculations  for  cohesive  materials.  Although  tensile  stresses  are 
consistent  with  the  assumed  conditions,  their  existence  in  limit-equilibrium 
slope-stability  solutions  may  lead  to  complications  and  errors.  Spencer's 
limit-equilibrium  procedure  (Spencer  1967)  encounters  various  degrees  of 
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numerical  instability  in  which  the  iterative,  trial  and  error  procedures  used 
to  calculate  the  factor  of  safety  may  not  produce  a  convergent  solution  when 
relatively  high  tensile  stresses  exist. 

30.  Tensile  stresses  resulting  from  high  cohesive-strength  components 
are  in  many  cases  unrealistic  and  may  easily  exceed  any  tensile-strength 
capacity  which  the  soil  can  provide  for  stability.  Often,  the  cohesive- 
strength  component  assigned  to  describe  the  strength  of  a  soil  simply  repre¬ 
sents  an  intercept  obtained  by  extrapolating  a  straight-line  failure  envelope 
to  zero  normal  stress  (a,  a)  on  a  Mohr-Coulomb  diagram.  While  such  a  cohesion 
intercept  implies  a  certain  tensile  strength,  tensile  strengths  are  seldom 
measured  or  considered  to  be  realistic  for  soils.  Accordingly,  tensile 
stresses  calculated  in  stability  analyses  employing  Mohr-Coulomb  parameters 
with  a  significant  cohesion  should,  in  many  instances,  be  considered  erroneous 
and  eliminated.  A  simple,  practical  solution  for  eliminating  tensile  stresses 
from  analyses  by  any  limit-equilibrium  procedure  is  to  introduce  a  vertical 
"crack"  and  to  terminate  the  upper  portion  of  the  shear  surface  at  a  vertical 
slice  boundary  with  an  appropriate  depth  below  the  crest  of  the  slope.  The 
depth  of  the  crack  d^  can  be  estimated  from  the  Rankine  active  earth- 
pressure  theory  as 


where 


(7) 


c  ,<p  =  mobilized  strength  parameters  (c  =  c/F  ,  tan  <p  -  [tan  <p]/F) 

mm  mm 

Y  =  unit  weight  of  soil 

Usually,  the  factor  of  safety  (c  and  p  )  can  be  estimated  with  a  reasonable 

m  m 

degree  of  accuracy  for  purposes  of  establishing  a  crack  depth  from  Equation  7 
prior  to  stability  calculations;  an  estimated  value  of  unity  for  F  is  often 
adequate.  Use  of  a  crack  depth  substantially  greater  than  the  depth  given  by 
Equation  7  should  be  avoided  because  it  will  not  only  eliminate  zones  of  ten¬ 
sile  stress  but  also  compressive  stress;  thus,  the  factor  of  safety  may  be 
overestimated.  Introduction  of  a  crack  with  an  appropriate  depth  eliminates 
both  unrealistic  tensile  stresses  and  virtually  all  numerical  instabilities 
which  are  associated  with  tensile  stresses  in  various  procedures. 
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Shear-Surface  Inclinations 


31.  A  second  potential  source  of  error  in  limit-equilibrium  analyses 
results  from  conditions  where  shear  surfaces  exit  steeply  upward  through  cohe¬ 
sionless  soil.  In  such  cases,  N  and  E  near  the  exit  portion  of  the  shear 
surface  may  become  very  large,  approaching  infinite  values  in  the  extreme,  or 
may  become  negative.  In  addition,  the  iterative  numerical  solution  used  in 
the  Simplified  Bishop  (Bishop  1955),  force-equilibrium,  and  complete  equilib¬ 
rium  procedures  to  calculate  the  factor  of  safety  may  oscillate  or  diverge. 
These  problems  are  often  associated  with  analyses  employing  circular  shear 
surfaces  for  embankments  on  foundations  where  a  relatively  thin  layer  of 
cohesionless  soil  is  underlain  by  a  thicker  stratum  of  weaker  soil,  often  soft 
clay.  In  these  cases,  the  critical  shear  surface  is  often  relatively  deep  and 
passes  steeply  upward  at  its  exit.  The  problem  of  very  large  or  negative 
stresses  near  the  toe  can  also  occur  for  relatively  shallow,  noncircular  shear 
surfaces,  which  have  a  relatively  flat,  horizontal  portion  at  a  shallow  depth 
exiting  abruptly  at  the  ground  surface.  When  the  assumed  inclination  of  the 
exiting  portion  of  the  surface  is  too  steep,  the  solution  for  the  factor  of 
safety  will  diverge,  or  unrealistic  values  for  the  factor  of  safety  and 
stresses  at  the  toe  of  the  shear  surface  will  result.  Whitman  and  Bailey 
(1967)  recognized  the  problem  associated  with  steeply  exiting,  circular  shear 
surfaces  and  discussed  the  problem  as  it  pertains  to  the  Simplified  Bishop 
procedure  (Bishop  1955)  although  the  problem  occurred  with  force-equilibrium 
and  complete  equilibrium  slice  procedures  as  well. 

32.  Wright  (1969)  has  examined  a  number  of  cases  where  unreasonably 
large  or  negative  values  have  been  calculated  for  the  stresses  near  the  toe  of 
the  shear  surfaces.  In  each  of  these  cases,  the  stresses  were  consistent  with 
those  calculated  by  trial  wedge  or  Coulomb  passive  earth-pressure  procedures. 
This  consistency  was  dependent  on  the  fact  that  the  shear-plane  and  interslice 
(earth-pressure)-force  inclinations  which  were  used  for  the  earth-pressure 
calculations  were  the  same  as  those  from  the  corresponding  limit-equilibrium 
slope  analyses.  However,  the  inclinations  of  the  shear  planes  which  were 
being  considered  did  not  represent  reasonable  critical  shear-plane  inclina¬ 
tions  and  would  not  have  been  expected  to  form  based  on  reasonable  passive 
earth-pressure  considerations;  instead,  the  inclinations  of  the  shear  planes 
were  much  steeper  than  passive  earth-pressure  theories  would  indicate. 
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33.  There  are  some  adjustments,  in  such  cases,  that  can  be  made  in  the 
side-force  assumptions  employed  in  the  various  limit-equilibrium  procedures. 
These  adjustments  can  improve  the  reasonableness  of  stresses  near  the  toe  of 
steeply  inclined  shear  surfaces.  However,  such  adjustments  may  require  tedi¬ 
ous  effort  on  the  part  of  the  user  and  may  only  partially  solve  the  problem. 

A  preferred  alternative  solution  is  to  adjust  the  shape  (inclination)  of  the 
shear  surface  in  the  cohesionless  material  to  conform  with  a  more  reasonable 
state  of  passive  shear  resistance.  An  appropriate  inclination  for  the  shear 
surface  can  be  estimated  from  Coulomb  passive  earth-pressure  theory.  The 
inclination  determined  in  this  manner  is  a  function  of  the  0  (earth-pressure 

resultant)  and  the  mobilized  friction  angle  $  of  the  cohesionless  soil. 

m 

The  shear-surface  inclinations  ol  are  plotted  versus  0  for  various  values 
of  <J>  in  Figure  6.  An  appropriate  a  can  be  estimated  from  this  figure  and 
the  shear  surface  adjusted  to  this  inclination.  Such  adjustments  to  a  will 
entirely  eliminate  unreasonable  stresses  in  the  passive  zone  and  provide  a 
distinctly  improved  solution. 


Figure  6.  Critical  ct's  from  Coulomb 
passive  earth-pressure  theory  for 
various  values  of  6  and  <f> 

(after  Wright  1 982 ) m 
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PART  IV:  CONCLUSIONS 


34.  In  conclusion,  the  relative  advantages  and  limitations  of  the  vari¬ 
ous  analysis  procedures  are  presented  in  Table  5.  These  procedures  differ 
principally  with  respect  to  the  static-equilibrium  conditions  which  they  sat¬ 
isfy  and  the  assumptions  which  are  made  to  achieve  a  statically  determinate 
solution  for  the  factor  of  safety.  The  Spencer  procedure  (Spencer  1967)  sat¬ 
isfies  all  static-equilibrium  requirements. 

Table  5 

Relative  Advantages  and  Limitations  of 
Limit-Equilibrium  Procedures 


Procedure 

Advantages 

Limitations 

Simplified  Bishop 

Simple;  value  for  F  is 
reasonably  correct 
calculations  can  be  per¬ 
formed  by  hand 

Restricted  to  circular 
shear  surfaces;  could 
give  incorrect  results 
for  pseudo-static 
analyses 

Force  equilibrium 
(Corps  of  Engineers 
Modified  Swedish 
side-force  assump¬ 
tion,  Lowe  and 
Karafiath's  side- 
force  assumption 
and  Corps  wedge 
side-force 
assumption 

Circular  and  noncircular 
shear  surfaces;  calcula¬ 
tions  can  be  performed  by 
hand 

Factor  of  safety  is  sen¬ 
sitive  to  side-force 
assumptions,  Lowe  and 
Karafiath's  side-force 
assumption  appears  to 
work  best;  solution 
does  not  satisfy 
moment  equilibrium 

Spencer 

Statically  correct;  circular 
and  noncircular  shear 
surfaces;  complete 
equilibrium 

Requires  use  of  computer 
for  most  practical 
applications 

35.  The  Simplified  Bishop  (Bishop  1955)  and  all  force-equilibrium  pro¬ 
cedures  do  not  fully  satisfy  the  requirements  of  static  equilibrium  and,  thus, 
are  not  statically  complete  procedures.  The  force-equilibrium  procedures  are 
relatively  inaccurate,  and  the  factors  of  safety  calculated  by  these  proce¬ 
dures  may  differ  by  50  percent  or  more  for  some  situations  from  the  calculated 
values  using  complete  equilibrium  procedures.  Accordingly,  the  usefulness  of 
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a  force-equilibrium  procedure  rests  primarily  in  its  suitability  for  hand  cal¬ 
culations  employing  noncircular  shear  surfaces.  The  Simplified  Bishop  proce¬ 
dure  is  probably  more  correct  than  either  of  the  force-equilibrium  procedures. 
Factors  of  safety  calculated  by  the  Simplified  Bishop  procedure  have  been 
found  to  agree  within  approximately  ±5  percent  with  those  calculated  by  com¬ 
plete  equilibrium  slice  procedures  for  relatively  simple  slope  and  loading 
conditions.  Thus,  for  at  least  some  cases,  the  Simplified  Bishop  procedure 
appears  to  be  relatively  correct.  However,  this  procedure  does  not  consider 
equilibrium  of  forces  in  the  horizontal  direction.  Accordingly,  use  of  the 
procedure  for  complex  slope  and  loading  conditions  where  significant  hori¬ 
zontal  external  loads  and  internal  body  forces  are  applied  may  lead  to  large 
inaccuracies . 

36.  Although  procedures  which  satisfy  all  requirements  for  static  equi¬ 
librium  are  statically  correct  procedures  and  most  appear  to  produce  essen¬ 
tially  the  same  value  for  the  factor  of  safety,  the  procedures  may  also 
produce  incorrect  or  unreasonable  values  for  the  factor  of  safety.  For  a 
variety  of  slopes,  the  geometry  assumed  for  the  shear  surface  in  the  "active" 
or  "passive"  zone  can  be  important  and  may  influence  the  computed  results. 
Introduction  of  a  vertical  crack  in  the  active  zone  or  adjustment  of  the 
inclination  of  the  shear  surface  in  the  passive  zone  is  often  required  to 
obtain  a  reasonable  solution.  The  assumption  of  a  constant  factor  of  safety 
along  the  shear  surface  implies  that  the  shear  strength  is  developed  simultan¬ 
eously  along  the  entire  shear  surface.  For  cases  where  this  is  not  true  due 
to  dissimilar  or  brittle  materials,  this  assumption  may  lead  to  incorrect  or 
unreasonable  results. 

37.  All  limit-equilibrium  analysis  procedures  are  subject  to  the  basic 
limitations  imposed  by  the  definition  of  the  factor  of  safety;  however,  all 
procedures  which  satisfy  complete  static  equilibrium  produce  almost  the  same 
value  for  the  factor  of  safety.  These  complete  equilibrium  procedures  provide 
a  useful  means  for  estimating  the  stability  of  earth  slopes  and  can  be  used 
with  considerable  confidence  for  their  intended  purposes. 
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APPENDIX  A:  LIMIT-EQUILIBRIUM  SLOPE-STABILITY  EQUATIONS  USED 
IN  THE  COMPUTER  PROGRAM  UTEXAS2  (Wright  1986)* 


Introduction 


1.  The  computer  program  UTEXAS2  has  the  capability  of  using  four  dif¬ 
ferent  limit-equilibrium  slope-stability  analysis  procedures  for  computing  the 
factor  of  safety  F  (Wright  1986).  The  four  procedures  are  Spencer’s  (1967) 
procedure;  Simplified  Bishop  (1955)  procedure;  force-equilibrium  procedure 
with  the  Corps  of  Engineers  Modified  Swedish  side-force  assumption,  Engi¬ 
neer  Manual  (EM)  1110-2-1902  (Headquarters,  Department  of  the  Army  1970);  and 
force-equilibrium  procedure  with  Lowe  and  Karafiath's  (1960)  side-force 
assumption.  Each  of  these  procedures  is  a  "procedure  of  slices"  in  which  the 
soil  mass,  bounded  by  the  surface  of  the  slope  and  an  assumed  shear  (sliding) 
surface,  is  divided  into  a  number  of  vertical  slices.  The  equilibrium  equa¬ 
tions  used  to  compute  the  factor  of  safety,  as  well  as  other  unknown  quanti¬ 
ties  computed  in  the  solution  by  each  of  these  limit-equilibrium  procedures, 
are  presented  in  this  appendix.  In  addition,  the  numerical  procedures  used  to 
solve  these  equations  in  the  computer  program  UTEXAS2  are  presented. 

2.  The  Modified  Swedish  and  Lowe  and  Karafiath  procedures  each  satisfy 
all  requirements  for  equilibrium  of  forces  on  individual  slices,  but  do  not 
satisfy  moment  equilibrium.  Except  for  differences  in  the  assumptions  made 
regarding  the  inclination  of  the  side  forces  in  these  two  procedures,  the  two 
procedures  are  otherwise  identical.  Accordingly,  presentation  of  the  equa¬ 
tions  and  numerical  solutions  is  combined  for  the  Modified  Swedish  and  Lowe 
and  Karafiath  (1960)  procedures  and  is  covered  under  the  single  heading  of 
"Force-Equilibrium  Procedures . " 

3.  The  general  nomenclature  and  sign  conventions  used  are  presented  in 
this  appendix  along  with  general  equations  and  derivations  common  to  all  pro¬ 
cedures.  The  remaining  paragraphs  of  this  appendix  cover  the  specific  equa¬ 
tions,  derivations,  and  numerical  solution  procedures  for  Spencer's,  the 
force  equilibrium,  and  Simplified  Bishop  procedures. 


*  References  cited  in  this  appendix  are  included  in  the  References  at  the  end 
of  the  main  text. 
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Nomenclature  and  Sign  Convention 

4.  A  typical  slice  and  the  system  of  forces  acting  on  the  slice  are 
shown  in  Figure  Al.  Coordinates  are  expressed  in  a  right-hand  Cartesian 
coordinate  system  with  the  x  axis  being  horizontal  and  directed  toward  the 
right,  and  the  y  axis  being  vertical  and  directed  upward  for  positive  values. 
Forces  are  considered  to  be  positive  when  they  act  in  the  positive  x  or  y 
direction,  unless  noted  otherwise.  Angles  are  measured  from  the  horizontal 
plane  (x  axis)  and  are  considered  to  be  positive  in  the  counterclockwise 
direction.  The  forces  acting  on  a  slice  consist  of  the  weight  of  the  slice 

W  ,  a  normal  and  a  shear  force  on  the  top  of  the  slice,  P  and  T  ,  respec¬ 
tively,  a  force  R  acting  on  the  base  of  the  slice  to  represent  internal 
reinforcement,  and  the  normal  and  shear  forces  on  the  base  of  the  slice  N 
and  S  ,  respectively.  Excluding  any  force  carried  by  reinforcement,  a  hori¬ 
zontal  force  KW  which  represents  the  body  force  for  seismic  loading  in  a 
"pseudo-static"  analysis,  and  the  normal  and  shear  forces,  E  and  X  , 
respectively  on  the  sides  of  the  slice,  these  forces  represent  the  forces  in 
the  soil.  The  normal  and  shear  forces  on  the  top  of  the  slice  represent  any 
loads  due  to  concentrated  forces  or  distributed  pressures  on  the  surface  of 
the  slope. 

5.  The  forces  shown  in  Figure  Al  represent  total  forces  including  any 

forces  due  to  water  pressures.  All  forces  are  shown  acting  in  the  direction 

in  which  they  are  assumed  to  act  when  their  value  is  positive.  Negative 

values  for  any  of  the  forces  act  in  the  opposite  directions  to  the  ones  shown 

in  Figure  Al.  The  weight  of  the  slice  is  assumed  to  act  through  the  midpoint 

of  the  slice,  and  the  normal  force  is  assumed  to  act  at  the  center  of  the  base 

of  the  slice  for  convenience.  Although  other  locations  might  be  assumed  for 

these  forces  (W  and  N) ,  the  locations  have  been  found  to  have  a  minor  effect 

on  the  factor  of  safety.  The  coordinates  of  the  center  of  the  base  of  the 

slice  are  designated  as  x  ,  y  ;  the  inclination  of  the  base  of  the  slice 

b  b 

is  by  the  angle  a  ;  the  location  of  the  normal  force  P  on  the  top  of  the 

slice  by  the  coordinates  x  ,  y  ;  the  inclination  of  the  top  of  the  slice 

P  P 

by  the  angle  8  ;  the  location  of  R  by  the  coordinates  x^_  ,  y^  on  the 
base  of  the  slice,  and  its  inclination  by  the  angle  ip  .  The  seismic  force 
acts  horizontally  along  a  line  with  the  y  coordinate  y^  . 
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6.  It  is  convenient  to  express  the  components  of  the  known  forces 

acting  in  the  horizontal  direction  on  a  slice  by  the  resultant  expressed 

as 

F,  =  -KW  +  P  sin  8  +  T  cos  6  +  R  cos  ip  (Al) 

n 

where  the  forces  acting  to  the  right  are  considered  positive.  Similarly,  the 
components  of  the  unknown  forces  acting  in  the  vertical  direction  are 
expressed  by  the  resultant  as 

Fv  =  -W  -  P  cos  8  +  T  sin  6  +  R  sin  ip  (A2) 

where  forces  acting  in  the  upward  direction  are  considered  positive.  The 
moment  produced  about  the  center  of  the  base  of  the  slice  by  the  known  forces 
is  expressed  by  as 

Mq  =  -P  sin  8(yp  -  yb)  -  P  cos  8(xp  -  xfe)  -  T  cos  8(yp  -  yb) 

+  T  sin  6 (xp  -  xfe)  +  KW(yfc  -  yfe)  -  R  cos  <Kyr  -  yb)  (A3) 

+  R  sin  i^(x  -  x,  ) 
r  b 

where  a  counterclockwise  moment  is  considered  to  be  positive. 

7.  The  shear  force  on  the  base  of  the  slice  can  be  expressed  as 

S  =  xA£  (A4) 

where 

t  =  average  shear  stress  on  the  base  of  the  slice 
AS,  =  length  of  the  base  of  the  slice 

The  average  shear  stress  is  related  to  the  shear  strength  and  factor  of  safety 
by  the  definition  of  the  factor  of  safety,  which  is  expressed  as 
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F  =  f  (A5) 

where  s  is  the  available  shear  strength. 

8.  The  shear  strength  can  be  expressed  in  terms  of  effective  stress 
shear-strength  parameters  c  and  $  or  total  stress  shear-strength  param¬ 
eters  c  and  <f>  .  In  the  case  where  effective  stresses  are  used,  the  avail¬ 
able  shear  strength  is  expressed  by  the  Mohr-Coulomb  equation  as 

s  =  c  +  (a  -  u)  tan  5  (A6) 

which,  when  substituted  into  Equation  A5  and  rearranged,  gives  the  following 
expression  for  the  shear  stress 


c  +  (g  -  u)  tan  <p 
F 


(A7) 


Finally,  combining  Equations  A7  and  A4,  the  following  expression  is  written 
for  the  shear  force  on  the  base  of  the  slice 


S  =  |  [cA£  +  (N  -  uA£)  tan  $] 

or,  since  A£  =  Ax  sec  a  ,  where  Ax  =  width  of  slice 

S  =  ^  [cAx  sec  a  +  (N  -  uAx  sec  a)  tan  <F] 

9.  In  the  case  where  the  shear  strength  is  expressed  using  total 
stresses, 


(A8) 


(A9) 


s  =  c  +  c  tan  $ 


(A10) 


Equations  A8  and  A9  become 


S 


=  —  [cA£  +  N  tan 


(All) 
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and 


S  =  4  [cAx  sec  a  +  N  tan  <j>] 

r 


(A12) 


respectively.  Comparison  of  Equations  All  and  A8  and  Equations  A12  and  A9 
shows  that  the  equations  for  total  stresses  may  be  arrived  at  from  the  equa¬ 
tions  based  on  effective  stresses  by  simply  setting  the  pore  water  pressure  to 
zero  and  replacing  c  and  <}>  in  the  effective  stress  equations  by  c  and 
<p  .  In  general,  all  of  the  equations  for  limit-equilibrium  slope-stability 
analyses  using  total  stresses  may  be  derived  from  the  corresponding  equations 
for  effective  stresses  by  setting  the  pore  water  pressure  u  to  zero  and 
replacing  c  and  <j>  in  the  effective  stress  equations  by  c  and  ((>  . 
Accordingly,  in  the  following  paragraphs  of  this  report,  all  equations  are 
presented  in  terms  of  effective  stress  shear  strength  parameters  with  the 
equations  for  total  stresses  being  an  implicit,  special  case,  i.e.,  u  ->•  0  ; 
c  -*  c  ;  4)  ->  <t>  . 


Spencer’s  Procedure 


10.  Spencer's  procedure  (Spencer  1967)  satisfies  all  requirements  of 
static  equilibrium.  The  side  forces  are  assumed  to  be  parallel,  i.e.,  they 
act  at  the  same  inclination.  The  value  of  the  inclination  0  is  considered 
to  be  an  unknown  which  is  solved  for  along  with  the  factor  of  safety. 


Derivation  for  the 
equilibrium  equations 


11.  In  deriving  the  equilibrium  equations  for  Spencer's  procedure,  it 
is  convenient  to  represent  the  side  forces  by  their  total  resultant  force  Q 
on  the  slice  as  shown  in  Figure  A2.  The  requirements  for  overall  equilibrium 
of  forces  in  the  horizontal  and  vertical  directions  are  then  expressed  by  the 
following  equations,  respectively: 


and 


EQ  cos  0=0 

(A13) 

EQ  sin  0=0 

(A14) 

where  the  summation  is  carried  out  for  all  slices.  However,  since  the  side- 
force  inclination  Q  is  assumed  to  be  constant  in  Spencer's  procedure. 
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Figure  A2.  Total  resultant  side  force  and 
line  of  action  coordinate  definition  for  a 
typical  slice 

Equations  A13  and  A14  are  reduced  to  the  following  single  equation  for  force 
equilibrium: 

IQ  -  0  (A15) 

12.  The  requirement  for  overall  equilibrium  of  moments  is  conveniently 
expressed  by  summing  moments  about  the  origin  of  the  coordinate  system  which 
gives 

IQ(x^  sin  0  -  y^  cos  0)  «  0  (A16) 

where,  again,  the  summation  is  performed  for  all  slices.  Expressions  for  Q 
and  y^  in  Equations  A15  and  A16  are  obtained  by  considering  the  equilibrium 
requirements  for  individual  slices. 

13.  Expression  for  Q  .  The  expression  for  Q  is  derived  by  first 
summing  forces  in  a  direction  perpendicular  to  the  base  of  the  slice.  Forces 
acting  in  an  upward  direction  toward  the  base  of  the  slice  are  considered 
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positive.  The  summation  for  forces  gives  the  following  equation  for 
equilibrium: 

N  +  F  cos  a  -  F  sin  a  -  Q  sin  (a  -  9)  =  0  (A17) 

v  h 

which  can  be  solved  for  the  normal  force  N  to  give 


cos  a  +  F,  sin  a  +  0  sin  (a  -  9) 
n 


(A18) 


Substitution  of  Equation  A18  into  Equation  A9  derived  for  the  shear  force  S 
from  the  Mohr-Coulomb  condition  gives  the  following  equation  for  the  shear 
force  on  the  base  of  the  slice: 


S  =  4  I  cAx  sec  a  +  [-F  cos  a  +  F,  sin  a  +  Q  sin  (a  -  9) 
r  |  v  h 

-  uAx  sec  a]  tan  <j>|  (A19) 

14.  An  expression  for  the  shear  force  on  the  base  of  the  slice  is  also 
obtained  by  resolving  forces  parallel  to  the  base  of  the  slice.  By  resolving 
forces  parallel  to  the  base  of  the  slice,  the  following  equilibrium  equation 
can  be  written: 


S  +  F  sin  a  +  F,  cos  a  +  Q  cos  (a  -  0)  =  0  (A20) 

v  n 

which  can  be  rearranged  to  give  the  following  expression  for  the  shear  force: 

S  =  -F  sin  a  -  F.  cos  a  -  0  cos  (a  -  9)  (A21) 

v  h 

The  expressions  for  the  shear  force  given  by  Equations  A19  and  A21  can  be 
equated  to  give 

4  IcAx  sec  a  +  f-F  cos  a  +  F  sin  a  +  Q  sin  (a  -  9) 

r  [  V  h 


-  uAx  sec  a]  tan  d>  =  -F  sin  a  -  F,  cos  a  -  Q  cos  (a  -  0)  (A22) 

v  n 
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Finally,  Equation  A22  can  be  solved  for  the  resultant  force  Q  to  give 


Q  = 


sin  a 


cos  a  -  ■=  Ax 
F 


sec  a  + 


(F  cos  a 
v 


F,  sin  a 
h 


+  uAx  sec  a) 


tan  <p 


ma 


(A23) 


where 


m  = 

a 


(A24) 


cos  (a  -  0)  +  sin  (a  -  6) 


tan  <p 


15.  Expression  for  y. 


The  line  of  action  of  the  resultant  force 


Q  is  expressed  by  the  coordinate  y^  located  on  the  line  of  action  at  a 
point  directly  above  the  center  of  the  base  of  the  slice.  The  coordinate 
y^  is  shown  in  Figure  A2.  Summing  moments  about  the  center  of  the 
base  of  the  slice  and  noting  that  the  sum  must  be  equal  to  zero  for  equilib¬ 
rium  gives 


-Q  cos  9(yq  -  yb)  +  Mq  »  0  (A25) 

where,  as  shown  earlier  by  Equation  A3,  represents  the  moment  about  a 

point  on  the  center  of  the  base  of  the  slice  due  to  all  of  the  known  forces 
(KW,  P,  R,  etc.).  Equation  A25  is  solved  for  y^  to  give 


M 

o 

^b  Q  cos  9 


(A26) 


Solution  of  equilibrium  equations 

16.  The  expressions  for  Q  (Equation  A23)  and  y^  (Equation  A26)  are 
substituted  into  the  equations  of  equilibrium  (Equations  A15  and  A16)  to  pro¬ 
duce  two  equations  in  two  unknowns  (F  and  0)  which  must  be  satisfied  to  sat¬ 
isfy  static  equilibrium.  The  solution  to  Equations  A15  and  A16  for  the  factor 
of  safety  and  side-force  inclination  is  accomplished  using  an  iterative  proce¬ 
dure  based  on  Newton's  method  for  solving  two  equations  in  two  unknowns. 
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For  assumed  values  of  the  factor  of  safety  Fg  and  side-force  inclination 
0q  ,  it  is  convenient  to  write  the  two  equilibrium  equations  in  the  form 

*  EQq  (A27) 

and 

R2  =  SQo(Xb  Sin  60  '  yQ0  C0S  6o)  (A28) 

where 

Qq  =  value  of  Q  based  on  the  assumed  values 

Fq.Sq.Rj,!^  =  force  and  moment  imbalances,  respectively,  based  on  the 
assumed  values  Fg  and  0g 

Application  of  Newton's  method  to  find  the  roots  to  Equations  A27  and  A28 
corresponding  to  R^  =  R^  =  0  ,  gives  the  following  for  the  new  estimates  for 
F  and  0  based  on  the  assumed  values: 

Fj  =  FQ  +  AF  (A29) 

and 


0j  -  e0  +  A0 


(A30) 


where  AF  and  A0  represent  adjustments  to  the  assumed  values  of  F  and 
0  ,  respectively,  to  be  used  for  the  next  iteration.  The  expressions  for  aF 
and  A0  are  as  follows: 


3R2 

3Rl 

R1 

30 

R2 

30 

3R2 

3R  j 

ro 

30 

3F 

3F 

30 

3Rj 

3R2 

R2 

3F 

'  R1 

W 

3Rj 

3R2 

3Rj 

3R2 

30 

3F 

3F 

30 

(A31) 


(A3  2) 


In  the  computer  program  UTEXAS2,  Equations  A31  and  A32  are  used  to  compute  the 
values  of  aF  and  A0  ,  respectively,  up  to  the  point  in  the  iterative  solu¬ 
tion  where  the  respective  values  become  less  than  0.5  and  0.15  radians.  Once 
the  values  become  less  than  these  limits  (aF  =  0.5,  A0  =  0.15  radians),  an 


A10 


"extended"  form  of  Newton's  method  is  used  based  on  the  following  two 
equations : 


3r 

R1  +  AF  3F 


3R1  l  2  j  2 

+  A03F-  +  2Af  7X  +  AFA6^e+iAe  zr 


3F 


38 


(A3  3) 


3R„ 

R2  +  AF  3F 


3R  1  2  3%  32r 

+  A9_  +  _Af  _2  +  AfA6^_ 


2  1  2  3  R2 
z  +  £  Ae^  — -=• 

2  302 


(A34) 


Equations  A33  and  A34  are  derived  from  Taylor  (1937)  series  expansions  includ¬ 
ing  the  second-order  terms.  Estimates  of  new  trial  values  are  obtained  by 
solving  these  two  equations  simultaneously  for  AF  and  A0  . 

17.  The  partial  derivatives  of  R^  and  R^  in  Equations  A31 
through  A34  are  obtained  from  Equations  A27  and  A28  and  are  as  follows: 


3ri  -  z  3Q 

If  2  If 


3Q 

10“  2  10 


32R. 


3f' 


-  2_ 

th 

3f36 


32R, 


30" 


=  E 


=  Z 


=  E 


2 

3_Q 

3F2 


2 

3_Q_ 

3F30 


2 

3  Q 
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(A35) 


(A36) 
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3R2 

W 


cos  6o) 


(A38) 


(A39) 


(A40) 


All 


sF  ■  £  H  xb(sin  6o  -  yQ0  cos  6o)  +  £Qo  (*b  cos  eo  +  y. 


Q0  8ln  90 


COS  6_  )  (A41) 


—T  =  2  “2  (Xb  Sin  e0  -  yq  COS  eo)  "  2E  If(  DF2  C0S  eO 

9F  9F  \  U  /  \  , 


-  ZQ0  l  —2  COS  90 

\9F  ) 


(A42) 


3  R2  32Q  / 

9F96  =  2  9F96  \Xb  sin  90  “  yQQ  COS  90 


+  2  If  ( 


3yO  \  90  /3yQ 

xb  COs  60  +  yQ  Sin  60  ‘  96"  C0S  60  /  '  2  96  \W  COS  6( 


/a2yQ  %  \ 

2Q0  \9F96  COS  e0  “  9F  sin  90/  (A43) 


92R  2 

— r  =  2  M  (xb sin  9o  -  yq cos  V 

96  96  v 


+  2Z  le  (xb  cos  eo  +  yq0  sin  90  "  96^  cos  90 


/  y0 
-  EO!*,  sin  6  -  y  cos  9n  -  2  ■—  sin  6  +  — —  cos  6 

\b  0  Q0  0  390  °  962 


o  ae. 


(A44) 


18.  In  evaluating  the  various  partial  derivatives  of  Q  in  Equa¬ 
tions  A35  through  A44,  it  is  convenient  to  write  the  expression  for  Q 
(Equation  A23)  as 


A12 


(A45) 


A13 


to\/D 


(A54) 


k  •  5>): 


where 


Mk  •*)(£•*£) 


2(c  +M(!S  +  1 9ca2' 
■ 2  \ci  +  r/  \9e"  +  ?  ir/  J 


tt  =  sir-  (a  ~ 6) 


(A55) 


=  -cos  (a  -  0)  tan 


(A56) 


=  -cos  (“  -  0) 


(A57) 


=»  -sin  (a  -  0)  tan  <p 


(A58) 


19.  Expressions  for  the  various  partial  derivatives  of  the  variable  y^. 
in  Equations  A40  through  A44  are  as  follows: 


(Q0  cos  0Q) 


I  M0  If  cos  eo 


(A59) 


90Q  *  /  \2  Mo(a0  COS  90  *  %  Sln  9o) 

\%  cos  e0' 


(A60) 


!!za  =  —±—  „  fe 

3F2  Q2  cos  0  °  10F2 


2  (^)2 

~  Q  \W 


(A61) 
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32y, 


3F30 


Q  _ 


-1 


Q0  cos  e0 


«  3Q  3Q  1_> 
3F  36  Qq j 


(A62) 


92y 


Q  _ 


36 


-1 


2  Mo 

Q0  cos  eQ 


2  ^  tan  6 
.  39 


0  “  g02  +  Q0  +  Q0  (30  “  Q0  tan  6o) 


(A63) 


Solution  for  remaining  unknowns 

20.  Once  thi  values  of  the  factor  of  safety  and  side-force  inclination 
are  determined  which  satisfy  the  equilibrium  Equations  A15  and  A16,  the 
remaining  unknowns  are  calculated.  The  remaining  unknowns  consist  of  the  nor¬ 
mal  force  on  the  base  of  the  slice  N  ,  the  side  force  Z  between  slices,  and 
the  locations  of  the  side  forces  y  .  The  normal  forces  are  calculated  from 
Equation  A18  which  was  derived  by  summing  forces  in  a  direction  perpendicular 
to  the  base  of  each  slice.  The  value  of  the  force  Q  in  Equation  A18  is 
calculated  from  Equation  A23.  Although  the  shear  force  is  not  actually  con¬ 
sidered  an  "unknown"  (it  is  known  if  N  and  F  are  known) ,  the  shear  force  can 
only  be  calculated  once  F  is  found;  the  shear  force  is  calculated  from 
Equation  A21. 

21.  Side  forces  Z  are  calculated  slice-by-slice,  beginning  with  the 

til 

first,  leftmost  slice.  In  general,  for  the  i  slice. 


where 


(A64) 


Z^  =  the  side  force  on  the  left  of  the  slice 
Z^+^  =  the  side  force  on  the  right  of  the  slice 

=  the  resultant  of  the  side  forces  on  each  side  of  the  slice 
The  side  forces  are  shown  in  Figure  A3.  The  resultant  of  the  side  forces 
is  calculated  for  each  slice  from  Equation  A23  once  values  for  the  factor  of 
safety  and  side-force  inclination  have  been  determined.  The  side  forces  are 
then  calculated  for  each  slice  beginning  with  the  first. slice.  For  the  first 
slice,  the  side  force  on  the  left  of  the  slice  Z^  must  be  zero  and,  thus. 


A15 


M 


(A65a) 

where 

Z2  =  the  force  on  the  right-hand  side  of  the  first  slice 
=  the  value  of  Q  for  the  first  slice 


Y 


1 — — - - - - - x 

Figure  A3.  Resultant  side  forces  and 
inclination  acting  on  the  sides  of  a 
typical  slice 

Application  of  Equation  A64  to  the  next  slice  then  gives 

Z3  =  Z2  “  Q2  (A65b) 

where 

Z3  =  the  side  force  on  the  right  of  the  second  slice 
Q,  =  the  value  of  Q  for  the  second  slice 
Equation  A64  is  applied  successively  to  the  remaining  slices  until  all  side 
forces  have  been  calculated. 

22.  The  location  of  the  side  forces  (line  of  thrust)  is  also  calculated 
slice-by-slice  beginning  with  the  first  slice.  By  summing  moments  about  the 
center  of  the  base  of  the  slice,  the  following  equation  can  be  written  for  any 
slice : 
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Mo  -  ZA  sin  0  -  Zi+1  sin  0  ^  -  Z±  cos  0(ytfi  -  yb> 

+  Zi+1  COS  0(yt>1+1  -  yb)  =  0  (A66) 

which  can  be  rearranged  to  give  the  following  expression  for  the  location  of 
the  side  force  on  the  right-hand  side  of  a  slice: 


yt,i+l  "  yb  - 


Mo  -  zi  stn  6  r  ~  zi+i  sin  6  r  "  zi  cos  e(yt,i  •  V 


Z  cos  0 
i+1 


(A67) 


For  the  first  slice  (i  =  1),  the  value  of  the  force  Z^  is  zero,  and  the 

value  of  the  location  y  .  is  of  no  significance  in  Equation  A67.  Thus, 

t  »  1 

y„  ...  ,  corresponding  to  the  location  of  the  side  force  on  the  right-hand 

side  of  the  first  slice,  can  be  calculated  once  Z^+^  has  been  calculated. 

This  process  can  be  repeated  for  the  next  slice  once  the  values  of  Z^  and 

Z  ,  and  the  value  of  y  .  have  been  calculated.  The  process  is  repeated 
I '  X  t )  1 

until  the  locations  of  all  of  the  side  forces  have  been  determined. 


Force-Equilibrium  Procedures 


23.  Two  force-equilibrium  procedures  are  used  in  the  computer  program, 
UTEXAS2.  The  first  of  these  procedures  uses  the  Corps  of  Engineers  Modified 
Swedish  side-force  assumption,  EM  1110-2-1902  (Headquarters,  Department  of  the 
Army  1970);  the  second  uses  the  Lowe  and  Karafiath’s  (1960)  side-force  assump¬ 
tion.  Both  procedures  satisfy  equilibrium  of  forces  in  the  vertical  and 
horizontal  direction,  but  do  not  satisfy  moment  equilibrium.  A  statically 
determinate  solution  is  obtained  by  assuming  the  inclination  of  the  resultant 
side  forces  between  slices.  In  the  case  of  the  Corps  of  Engineers  Modified 
Swedish  side-force  assumption,  the  side  forces  are  assumed  to  be  parallel  (all 
side  forces  have  the  same  inclination),  and  the  inclination  must  be  selected 
and  input  by  the  user.  According  to  EM  1110-2-1902  (Headquarters,  Department 
of  the  Army  1970),  the  side-force  inclination  would  normally  be  taken  to  be 
equal  to  the  "average  embankment  slope"  although  in  actual  practice  other 
inclinations  might  be  assumed.  In  the  case  of  Lowe  and  Karafiath's  side-force 
assumption,  the  side  forces  are  assumed  to  act  at  the  average  of  the 
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inclination  of  the  slope  (or  ground)  surface,  directly  above,  and  the  shear 
surface,  directly  below,  each  vertical  boundary  between  slices.  The  side- 
force  inclinations  vary  from  slice  to  slice. 

Derivation  of  the  equilibrium  equations 

24.  In  deriving  the  equations  used  to  compute  the  factor  of  safety  by 
the  force-equilibrium  procedures,  it  is  convenient  to  express  the  forces  on 
each  side  of  the  slice  by  the  resultant  forces  Z^  and  Z^+^  acting  on  the 
left  and  right  side  of  the  slice,  respectively,  and  the  respective  inclina¬ 
tions  0  and  0  .  The  resultant  forces  are  illustrated  in  Figure  A3. 

Summation  of  forces  in  the  vertical  direction  gives  the  following  equilibrium 
equation: 


F  +  Z,  sin  0.  -  Z.j.  sin  0  .  +  N  cos  a  +  S  sin  a  =  0  (A68) 

v  i  l  l+l  l+l 

Similarly,  the  summation  of  forces  in  the  horizontal  direction  produces  the 
following  equilibrium  equation: 

F,  +  ?■  cos  0.  -  Z  ,  cos  8.  .  -  N  sin  a  +  S  cos  a  =  0  (A69) 

hi  l  i+1  i+1 

Substituting  the  expression  for  the  shear  force  S  given  by  Equation  A68  into 
Equations  A68  and  A69,  respectively,  gives  the  following  two  equations: 


F  +  Z,  sin  0 
v  i  i 


Zi+1  Sin  6i+l  + 


and 


F,  +  Z,  cos  0 
hi  i 


Zi+1  COS  8i+l  + 


N^cos  a  +  ta^— £  sin 

+  ^c  -  u  tan  sin  a  =  0 

N^sin  a  +  t3p~  ^  cos  aj 

+  ^c  -  u  tan  ?  j  ^  cos  o  =  0 


(A70) 


(A71) 


Equations  A70  and  A71  can  be  combined  to  eliminate  the  unknown  normal  force 
N  and  solved  for  the  side  force  acting  on  the  right  side  of  the  slice 

to  give 
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cos  (a  -  0i+1)  +  --a-p  -  sin  (a  -  @i+1) 


F  sin  a  +  F,  cos  a 
v  h 


+  Z.  cos  (a  -  0  )  +  — aa  --  [-F  cos  a  +  F,  sin  a  -  uAJt 
i  i  F  v  h 


+  sin  (a  -  0^) ]  + 


(A72) 


Solution  of  equilibrium 
equations  for  the  factor  of  safety 


25.  The  solution  for  the  factor  of  safety  is  obtained  using  an  itera¬ 

tive  procedure  based  on  Newton's  method.  A  factor  of  safety  is  assumed,  and 
the  side  forces  are  computed  slice-by-slice  using  Equation  A73.  Beginning 
with  the  first  slice  and  noting  that  for  the  first  slice  is  zero,  the 

force  Z.+^  on  the  right  of  the  slice  is  calculated.  The  force  z±+i  calcu¬ 
lated  for  the  right  of  the  first  slice  becomes  the  force  Z^  on  the  left  of 
the  next  slice,  and  the  force  on  the  right  of  the  next  slice  can  then  be 
calculated  once  again  using  Equation  A72.  This  procedure  is  repeated 
slice-by-slice  to  the  last  slice  where  the  force  on  the  right  side  of  the  last 
slice  is  calculated.  If  the  force  on  the  right  of  the  slice  is  acceptably 
small,  the  factor  of  safety  is  considered  to  be  correct,  and  the  remaining 
unknowns  can  be  calculated  as  described  in  the  following  paragraphs.  Other¬ 
wise,  a  new  value  is  assumed  for  the  factor  of  safety,  and  the  process  is 
repeated  until  the  force  on  the  right  of  the  last  slice  is  acceptably  small. 

26.  In  the  computer  program  UTEXAS2,  an  iterative  procedure  based  on 
Newton's  method  is  used  to  compute  the  factor  of  safety.  A  factor  of  safety 
Fq  is  assumed,  and  a  new  estimate  for  the  factor  of  safety  F^  is  obtained 


F1  =  F0  +  AF 


(A73) 


where 


(A74) 
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The  values  of  and  3Zi+1/3F  in  Equation  A74  are  evaluated  for  the 

last,  rightmost  slice  using  the  assumed  value  for  the  factor  of  safety  . 

In  evaluating  the  partial  derivative  3Z^+j/3F  ,  it  is  convenient  to  write 
Equation  A72  in  the  form 


Z 


i+1 


(A75) 


where 


C,  =  F  sin  a  +  F,  cos  ct  -f  Z .  cos  (<*  -  0  ) 
1  v  h  1  i 


C.  =  tan  $  [-F  cos  a  +  F,  sin  -  uAi  +  Z  sin  (a  -  9,)]  + 
2  v  h  i  i 


cAi 


=  cos(«  -  e1+1) 


=  tan  ?  sin  (<*  -  0i+1> 


The  partial  derivative  3Zi+^/3F  then  becomes 


3Z 


i+1 


3F 


=  < 


C3  +F 


3ci  1  /  3C2  ^ 

3F  +  p2  (F  3F  '  C2, 


where 


C2  \  C4 

+  ,ci+r)7 


-2 


>  C3  +  r 


3ci  3zi 

Tr  =  7T  cos(a  '  V 


3c2  3zi 

W~  m  W  sin(a  "  V  tan  * 


(A76) 

(A77) 

(A78) 

(A79) 


(A80) 


(A81) 


(A82) 


Equation  A6  is  used  to  compute  successive  trial  values  for  the  factor  of 
safety  until  the  changes  on  successive  trials  Af  and  the  force  imbalance 
Z  on  the  last  slice  become  acceptably  small. 
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Solution  for  remaining  unknowns 

27.  Once  the  value  of  the  factor  of  safety  is  determined  by  satisfying 

the  equilibrium  requirement  that  the  side  force  must  be  essentially 

zero  for  the  last  (right-most)  slice,  the  remaining  unknowns  are  calculated. 
The  remaining  unknowns  consist  of  N  on  the  base  of  the  slice  and  Z  between 
slices.  However,  the  side  forces  are  calculated  as  part  of  the  iterative  pro¬ 
cedure  used  to  calculate  the  factor  of  safety,  and,  thus,  only  N  on  the  base 
of  each  slice  remains  to  be  calculated.  The  expression  used  to  calculate  the 
normal  forces  is  obtained  by  resolving  forces  in  a  direction  perpendicular  to 
the  base  of  the  slice  and  by  solving  the  resulting  equilibrium  equation  for 
the  normal  force  to  give 

N  =  -F  cos  a  +  F,  sin  a  +  Z.  sin(a  -  0.)  -  Z, sin(a  -  0  )  (A83) 

v  h  i  i  i+1  i+r 

Although  the  shear  force  is  not  actually  considered  an  unknown  (it  is  known  if 
N  and  F  are  known) ,  the  shear  force  can  only  be  calculated  once  F  is 
found;  the  expression  used  to  calculate  the  shear  force  is  derived  by  summing 
forces  in  a  direction  parallel  to  the  base  of  the  slice  and  by  solving  the 
resulting  equilibrium  equation  for  the  shear  force  to  give 

S  =  -F  sin  a  -  F  cos  a  -  Z,  cos(a  -  0  )  +  Z..,  cos(a  -  0  )  (A84) 

v  n  i  l  l-t-i  x  •  I 

Simplified  Bishop  Procedure 

28.  The  Simplified  Bishop  procedure  (Bishop  1955)  is  based  on  the 

assumption  of  a  circular  shear  surface.  Side  forces  are  assumed  to  act  in  the 
horizontal  direction;  i.e.,  there  is  assumed  to  be  no  shear  force  between 
slices.  The  Simplified  Bishop  procedure  satisfies  equilibrium  of  forces  in 
the  vertical  direction  for  each  slice  and  equilibrium  of  moments  about  the 
center  of  the  circular  shear  surface  for  the  entire  free  body  composed  of  all 

slices  (overall  moment  equilibrium).  This  procedure  has  been  extended  by 

Wright  (1986)  to  account  for  external  loads. 

Derivation  of  the  equilibrium  equations 

29.  An  expression  for  the  normal  force  on  the  base  of  each  slice  is 
obtained  first  by  summing  forces  in  the  vertical  direction.  The  resulting 
equation  for  equilibrium  of  forces  in  the  vertical  direction  is 
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Solution  of  equilibrium 
equation  for  the  factor  of  safety 

31.  Newton's  method  is  used  in  the  computer  program  UTEXAS2  to  solve 
Equation  A91  for  the  factor  of  safety.  For  an  assumed  factor  of  safety  FQ  , 
Equation  A91  is  written  as 


Mi  =  lMo  +  ZFv(xb  -  xc)  -  EFh(yb  “  yc} 


+  —  £[cAx  +  (-F  -  uAx)  tan  ?]  m  (A92) 

F„  v  a 


where  represents  the  moment  imbalance  based  on  the  assumed  factor  of 

safety.  The  new  estimate  for  the  factor  of  safety  is  written  as 


Fj  -  F0  +  4F 


(A93) 


where,  by  Newton's  method,  &F  is  expressed  as 


(A94) 


In  computing  the  partial  derivative  9m/9F  ,  it  is  convenient  to  write 
Equation  A92  as 

C1 

Mi  =  EMo  +  EFv(xb  -  xc}  "  SFh(yb  -  yc}  +  RZ  F0C2 

where 


(A95) 


C,  =  c^x  -  (F  +  u&x)  tan  <f> 
1  v 


(A96) 


C2  =  cos  o 


(A97) 


=  sin  o  tan  ? 


(A98) 
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Then 


3M 

aF 


=  -RI 


<C2Fo  +  C3) 


2  C2 


(A99) 


Solution  for  remaining  unknowns 

32.  The  only  unknown,  in  addition  to  the  factor  of  safety,  which  is 
calculated  in  the  Simplified  Bishop  procedure  (Bishop  1955)  is  the  normal 
force  on  the  base  of  the  slice.  The  expression  used  to  calculate  the  normal 
force  is  obtained  by  substituting  the  expression  for  the  shear  force  (Equa¬ 
tion  A9) ,  obtained  from  the  Mohr-Coulomb  equation,  into  the  expression  for  the 
normal  force  (Equation  A86) ,  obtained  by  resolving  forces  in  the  vertical 
direction,  which  gives 

N  =  |  -F  -  fcAx  sec  a  +  (N  -  uAx  sec  a)  tan  <ji]  sin  a>  (A100) 

COS  Ot  I  r  I 

Rearranging  Equation  A100  gives  the  following  equation  for  the  normal  force  on 
the  base  of  each  slice: 

N  =  |-FV  -  |r  (c  -  u  tan  $)Ax  tan  a  m^J  (A101) 

where  m^  is  defined  in  Equation  A89. 


APPENDIX  B:  NOTATION 

Total  stress  equi/alent 

Shear-strength  parameter  in  terms  of  effective  stresses 
Mobilized  strength  parameter 
Depth  of  crack 

Horizontal  side  forces  between  slices 
Horizontal  forces  on  right  boundaries  of  slices 
Horizontal  forces  on  left  boundaries  of  slices 
Factor  of  safety 

Assumed  values  of  factor  of  safety 

Resultant  expressing  components  of  known  forces  acting  in  the 
horizontal  direction  on  a  slice 

Resultant  expressing  components  of  known  forces  acting  in  vertical 
direction  on  a  slice 

New  estimate  for  the  factor  of  safety 

Location  of  side  forces  on  the  side  of  the  slice 

Water  forces  on  side 

Water  forces  on  left  side 

Water  forces  on  right  side 

Horizontal  force  representing  body  force  for  seismic  loading  in  a 
pseudo-static  analysis 

Location  of  normal  force  on  base  of  slice 

Moment  produced  above  center  of  base  of  slice  by  known  forces 

Number,  of  slices 

Normal  forces  on  base  of  slice 

External  surface  loads 

Total  resultant  force  of  side  forces 

Value  of  Q  for  first  and  second  slices 

Bi 


Force  acting  on  base  of  slice  to  represent  internal  reinforcement 

Available  shear  strength 

Shear  force  on  base  of  slice 

External  shear  forces 

Pore  water  pressure 

Water  forces  on  base  of  slice 

Weight  of  the  slice 

Coordinate  of  the  center  of  base  of  slice 

Shear  forces  on  right  boundary  of  slice 

Coordinate  on  top  of  slice 

Coordinate  on  base  of  slice 

Vertical  shear  forces  on  slice  boundary 

Vertical  shear  forces  on  left  boundary  of  slice 

Coordinate  of  the  center  of  base  of  slice 

Cooridnate  of  the  centroid  of  the  slice 

Coordinate  on  top  of  slice 

Coordinate  on  base  of  slice 

Cooridnate  of  side  force  location 

Coordinate  of  resultant  force 

Resultant  side  force 

Resultant  side  force  on  left  of  slice 

Resultant  side  force  on  right-hand  side  of  the  first,  second,  and 
third  slices 

Shear-surface  inclinations 
Ground-surface  inclination 
Unit  weight  of  soil 
Adjustment  to  assumed  value  of  F 


Length  of  slice  base 


Ae 

e 

eo 

a 


V 

♦ 

^m 

T 


Adjustment  to  assumed  value  of  e 

Inclimation  of  resultant  interslice  force 

Assumed  value  of  side-force  inclination 

Total  normal  stress 

Effective  normal  stress 

Total  stress  equivalent 

Modilized  strength  parameter 

Shear-strength  parameter  in  terms  of  effective  stresses 
Mobilized  friction  angle  of  cohesionless  soil 
Shear  stress  required  for  just-stable,  static  equilibrium 
Inclination  of  reinforcement  force 
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Report  2  Simplified  Frame  Model  (Phase  D) 

Report  3:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies— Open  Section 

Report  4:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies— Closed  Sections 


Date 

Jan  1983 
Jun  1983 
Jul  1983 

Sep  1983 
Sep  1983 

Oct  1983 
Jan  1984 

Aug  1984 

Sep  1984 

Sep  1984 

Oct  1984 

Jun  1986 

Jan  1987 
Apr  1987 
May  1987 
May  1987 
Jun  1987 

Jun  1987 
Aug  1987 


(Continued) 


WATERWAYS  EXPERIMENT  STATION  REPORTS 
PUBLISHED  UNDER  THE  COMPUTER-AIDED 
STRUCTURAL  ENGINEERING  (CASE)  PROJECT 

(Concluded) 
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Instruction  Report  GL-87-1 

Instruction  Report  ITL-87-5 
Instruction  Report  ITL-87-6 

Technical  Report  ITL-87-8 
Instruction  Report  ITL-88-1 
Technical  Report  ITL-88-1 
Technical  Report  ITL-88-2 
Instruction  Report  ITL-88-2 
Instruction  Report  ITL-88-4 
Instruction  Report  GL-87-1 


Title 

Finite  Element  Studies  of  a  Horizontally  Framed  Miter  Gate 
Report  5:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies— Additional  Closed  Sections 
Report  6:  Elastic  Buckling  of  Girders  in  Horizontally  Framed 
Miter  Gates 

Report  7;  Application  and  Summary 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  I. 

User's  Manual 

Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Criteria  Specifications  fQr  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 
Gates  (CMITER) 

Procedure  for  Static  Analysis  of  Gravity  Dams  Using  the  Finite 
Element  Method  —  Phase  la 

User’s  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

User's  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

User's  Guide  for  Design  and  Investigation  of  Horizontally  Framed 
Miter  Gates  (CMITER) 

User's  Guide  for  Revised  Computer  Program  to  Calculate  Shear. 
Moment,  and  Thrust  (CSMT) 

User's  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  II, 
Theory 


Oate 

Aug  1987 

Aug  1987 

Oct  1987 
Dec  1987 

Jan  1988 
Feb  1988 
Apr  1988 
Apr  1988 
Jun  1988 
Sep  1988 
Feb  1989 


